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/  An  analytical  study  using  ths  STAGSC-1  coaputsr  cods 
Mas  conductsd  to  dstsrsina  ths  effects  of  notchss  and 
unsynsstric  boundary  conditions  on  ths  load  bsaring 
capability  and  radial  displacsasnts  of  axially  loadsd 
cylindrical  pansls.  Ths  graph its— spoxy  panel  consisted  of 
an  8  ply  laminate  wi th  symmetric  ply  orientations.  A  finite 
element  mesh  refinement  in  the  vicinity  of  the  cutout  mas 
conducted  in  order  to  study  the  non-linear  collapse 
analysis.  In  additon,  the  effect  of  different  size  cutouts 
Mith  aspect  ratios  <axial  dimension  of  cutout  divided  by 
circumf erential  dimension)  of  2.0  or  0.5  Mas  studied. 
Finally,  analytic  results  for  a  2"x2"  cutout  Mith  a  varying 
set  of  boundary  conditions  Mere  compared  to  experimental 
findings. 

It  Mas  found  that  as  the  surface  area  of  the  cutout 
increased,  the  buckling  load  decreased.  Cutouts  Mith  aspect 
ratios  of  2.0  Mere  found  to  be  capable  of  carrying  higher 
loads  than  notches  Mith  an  aspect  ratio  of  0.5.  By 
comparing  analytic  results  Mith  experimental  data,  it  Mas 
found  that  unsymmetric  boundary  conditions  for  the  2"x2" 
cutout  better  approximated  experimental  data  than  Mhen 

symmetric  boundary  conditions  Mere  considered. ^ . It  Mas  also 

determined  that  some  of  the  discrepancy  between  experimental 
and  analytic  results  can  be  attributed  to  initial 
eccentricities  introduced  into  the  panel  by  the  test 
apparatus. 


Background 


Tht  strength  and  stability  of  aircraft  and  space 
vehicles  are  influenced  by  tha  prasanca  of  cutouts  or 
notchas  in  thair  skins.  Many  vehiclas  in  tha  NASA  spaca 
prograa  ara  charactarizad  by  large,  lightwaight  shall 
structuras  that  can  ba  aithar  stringer  stiffened  or 
unstiff an ad,  and  can  hava  aithar  rainf oread  or  n on- 
reinforced  cutouts.  These  shall  structuras  can  ba  analyzed 
by  finite  eleaent  techniques.  Anderson  tl3  discusses  hoe 
finite  element  techniques  hava  been  used  in  tha  Apollo 
launch  vehicle  analysis  and  Skylab  stability  analysis. 
Design  studies  have  bean  conducted  using  finite  element 
techniques  on  the  Spaca  Shuttle  oxygen  tanks  and  cargo  doors 
C23.  Skogh  C31  describes  the  use  of  finite  element  analysis 
in  the  Skylab  structural  design. 

Of  particular  interest  in  space  vehicle  design  is  how 
structural  stability  is  affected  by  cutouts  in  the  shell 
structure.  Very  few  references  can  be  found,  however, 
dealing  with  axially  loaded  composite  cylindrical  panels 
with  cutouts.  Considerable  work  has  been  done  on  isotropic 
plates  and  cylindrical  panels  E4-73.  Sob el  C63  and  Almroth 
Ml  also  studied  how  varying  sets  of  boundary  conditions 
affected  the  buckling  behavior  of  unstiffened  panels.  Sobel 
C63  found  that  preventing  circumf erential  edge  displacement 
is  the  most  important  in-plane  boundary  condition  from  the 


point  of  view  of  increasing  buckling  load.  Several  authors 

133  hava  invastigatad  the  buckling  behavior  and  stress 
states  of  isotropic  panels,  with  cutouts,  under  axially 
coepressive  loads.  These  studies  were  done  nueerically 
generally  using  finite  difference  techniques.  Plates  and 
panels  with  reinforced  cutouts  have  also  been  analyzed  C133. 
Most  tests  have  been  conducted  using  circular  cutouts,  but 
sooie  investigations  have  been  done  with  elliptic  and  square 
cutouts. 

In  recent  years,  there  have  been  numerous  technological 
advances  in  the  application  of  fiber— reinforced  coeposites 
to  aircraft  and  space  vehicles.  Research  into  composite 
cylindrical  panels  without  cutouts  has  been  carried  out  in 
the  past  C 14-223,  but  relatively  little  work  can  be  found 
dealing  with  composite  panels  with  cutouts.  Satyamurthy 
C203  found  that  boundary  conditions,  load  systems,  location 
of  concentrated  loads,  and  material  properties  influence  the 
nature  of  the  buckling  behavior  in  composite  panels.  These 
types  of  analyses  must  also  be  applied  to  composite 
cylindrical  panels  with  cutouts  in  order  to  more  fully 
understand  the  effect  that  a  cutout  has  on  buckling 
behavior. 

Janisse  1233  is  one  of  the  few  authors  found  to  address 
the  problem  of  cutouts  in  composite  cylindrical  panels. 
Singer  C243  has  consolidated  much  of  the  work  done  on 
buckling  of  shells  into  a  concise  report.  This  report 
encompasses  isotropic  shells,  stiffened  shells,  and 
composite  shells  and  addresses  special  problems  and 
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techniques  used  in  shell  analysis. 

Composite  panels  have  been  studied  in  greater  detail 
froe  an  analytic  point  of  view  than  they  have  been 
experieentally.  This  is  due  to  eanuf acturing  difficulties 
and  lack  of  experieental  test  fixtures  needed  to  apply  a 
purely  coepressive  load.  Analytic  problees  are  eore  coeplex 
when  working  with  coeposite  Materials  because  of  their 
anisotropic  Material  properties  and  the  non  hoeogeneous 
coaplex  fiber /Matrix  interactions.  Nuaerous  coaputer  codes 
have  been  developed  to  analyze  these  types  of  structures. 
The  Structural  Analysis  of  General  Shells  CSTABSC— 1) 
coaputer  prograa  was  used  for  finite  eleaent  analysis  in 
this  thesis. 

Since  coaposites  are  used  alaost  exclusively  in 
aerospace  structures,  ao re  work  needs  to  be  done  in  the  area 
of  cylindrical  panels  with  cutouts  in  order  to  More  fully 
understand  the  effect  of  the  cutout  and  how  varying  the 
boundary  conditions  affects  the  buckling  behavior  of  a 
panel.  This  thesis,  which  is  a  continuation  of  the  work 
done  by  Janisse  C231,  addresses  the  issue  of  saall  cutouts 
in  composite  cylindrical  panels  and  the  effect  of 
unsyamctric  boundary  conditions  on  buckling  behavior. 

PfrjpeMVff 

The  purpose  of  this  thesis  is  threefold.  First,  the 
effect  that  cutouts  or  notches  of  various  aspect  ratios  have 
on  the  buckling  behavior  of  composite  panels  will  be  studied 


analytically  using  the  STA6SC-1  coaputar  p^ograa.  For  tha 
purposes  of  this  thesis,  a  notch  will  ba  dafinad  as  a  cutout 
of  lass  than  51  of  tha  total  panal  araa.  Second,  tha 
affacts  of  unsyaaatric  boundary  conditions  on  a  panal  with  a 
squara  cutout  will  ba  anal yz ad.  Finally,  analytic  rasults 
will  ba  compared  to  axpariaantal  data.  This  thasis  will 
anabla  a  battar  overall  understanding  of  tha  nature  and 
stability  of  coaposite  cylindrical  panels  under  axial  loads. 

Approach 

Tha  STA6SC-1  finite  alaaant  coaputar  prograa  C263 
developed  at  Lockheed  was  used  to  study  the  collapse  of 
cylindrical  panels  with  cutouts.  The  analysis  carried  out 
is  nonlinear  due  to  the  presence  of  a  cutout.  Tha  nonlinear 
collapse  analysis  branch  in  STAGSC-1  uses  large  displacaaant 
theory  and  Moderate  rotations  along  with  a  eodified  Newton— 
Raphson  iteration  scheae. 

Six  different  size  cutouts  located  at  tha  canter  of  tha 
panal  were  investigated.  A  panal  with  a  2"x2*  notch  (a 
cutout  coaprising  lass  than  51  of  panal  area)  was  used  as  a 
basis  for  a  aash  refinement  study.  These  results  were  then 
applied  to  the  finite  eleaent  aesh  design  for  the  other  five 
cutouts.  Of  these  six  cutouts,  two  were  square,  coaprising 
2.7%  and  11%  of  the  surface  araa  of  tha  panel,  respectively. 
Tha  other  four  cutouts  were  rectangular  with  aspect  ratios 
of  either  0.5  or  2.0.  Aspect  ratio  is  defined  as  the  axial 


diaension  of  the  cutout  divided  by  the  circumferential 


11. THEORY  AND  HOPFi  i  tub 


Clmltll  Laai nation  Theory 


In  order  to  understand  coaposite  laeinae  behavior, 
basic  constitutive  relationships  for  an  individual  laeina 
oust  be  explained.  The  reader  should  refer  to  References 
C143  and  E263  for  an  in-depth  developeent  of  these  relations 


The  stress-strain  expressions  for  an  individual  laeina 
are  essentially  the  stress— strain  relations  for  plane  stress 
in  an  orthotropic  aaterial.  The  fol lowing  relations  are  the 
basis  for  the  stiffness  and  stress  analysis  of  an  individual 
laeina  subjected  to  farces  in  its  own  plane.  For  a  laeinae 
of  orthotrapic  aaterial,  the  stress-strain  equations  in 
principal  aaterial  coordinates  are 


Q  q  e 

11  12 
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where  and  are  noreal  strains  and  yi2  is  the  shear 

strain  for  the  laeina  principal  axes  as  in  Figure  2.1,  and 
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Figure  2.1.  Definition  of  Coordinate  Systees 


where  and  E2  are  Young's  eoduli  in  the  1,  and  2 

directions,  respectively.  is  Poisson's  ratio  for 

transverse  strain  in  the  j  direction  when  loaded  in  the  i 
direction  and  G^2  is  the  shear  eodulus  in  the  1-2  plane. 
Therefore, 


(3) 


The  preceding  stresses  and  strains  Mere  defined  in  the 
principal  Material  directions  for  an  orthotropic  eaterial. 
If  the  principal  directions  do  not  coincide  with  coordinate 
directions,  stress  and  strain  transformations  are  used  to 
define  the  following  stress-strain  relations  in  xy 


coordinate 
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Figure  2.2. 
Therefore, 


'try  of  an  N  Layered  Laeinate 
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Mhere  the  eaterial  coef f icientc,  QijP  are  eyeeetric  and 
defined  as 

Q11  »  Q^co*4#  ♦  2CQ12+2Q66>sin^9cos2e  +  Q^Bin4^ 

012  -  (01 1+022-4066) si n20cos20  ♦  Q12(sin40+cos40  > 

6^  -  QjjSin4©  ♦  2(O12+2Q66>sin20cos20  +  Q22cos40  (5) 

61&  -  CQ11-Q12-'2O&6>sin0cos30  +  CQl2”O22+2Q6&>sin30cos0 
Q26  -  (Q1i~Oi2“2Q66>*in30co*0  ♦  CQ12-Q22+2Q&6)sin0cos30 
0^  -  <0l  1+Q22_2Q12~2Q66>  *in20cos20  ♦  <sin40+co*40  > 

For  a  detailed  derivation  of  the  transforeed  reduced 
stiffnesses,  see  Reference  €263. 

In  order  to  define  the  stress  and  strain  variations 
through  the  thickness  of  a  laminate,  stress-strain  behavior 
of  an  individual  lamina  Hill  be  expanded.  For  the  k**1  layer 


of  a  Multilayered  lMinatt,  Equation  <4>  can  be  Mr  it  tan  as 


C"k  “  C«k  COk  <*> 

Refer  to  Figure  2.2  for  tha  gaoaatry  of  an  N  layarad 
laainata  of  thicknass  t  (  aiddla  surfaca  location  and 
location  of  tha  kth  laaina. 

It  is  assuaad  that  aach  lay ar  can  ba  characterized  as 
an  orthotropic  Material  Mhosa  properties  are  known  with 
respect  to  its  Material  axes.  Tha  laainata  is  assuaad  to 
consist  of  perfectly  bonded  laainaa.  Tha  displaceaents  are 
continuous  across  laaina  boundaries  so  that  no  laaina  can 
slip  relative  to  one  another .  Also,  for  a  thin  laainata, 
noraal  sections  reaain  n areal  to  tha  aiddla  surface  after 
bending.  These  assuaptiona  for  a  laainata  constitute  tha 
Kirchoff-Love  hypothesis  for  shells  [263. 

It  is  necessary  to  develop  the  nonlinear  displaceaent 
cylindrical  shall  theory  in  order  that  the  reader  obtain  an 
appreciation  of  the  overall  theoretical  approach.  The 
aiddla  surface  strain— displaceaent  relations  for  a  laainata 
with  Moderately  large  rotations  of  tangents  to  the  panel 
reference  surface  are  given  by 

€k°  -  u°,x  -  u,K  ♦  1/2  4>2  *  1/2  <f>2 

€y°  -  v°,  -  v,y  «■  W/R  ♦  1/2  4>2  -  1/2  <t>2  (7) 

yxy°  m  v’x  +  u,y  ♦ 

ah  ere  u,  v  and  w  denote  the  axial,  circuaferential  and 

radial  components  of  the  aidsurface  displacement  and  4>  , 

x  y 

and  ^  are  the  components  of  rotation  about  the  coordinate 
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Figure  2.3.  Positive  Farces  end  Moments  Acting  on  a  Laminate 


lines  and  about  the  noreal  to  the  surface.  The  eiddle 
surface  curvatures  are 
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In  teres  of  displaceaents»  the  rotation  components  are 


0y  -  -W»y  ♦  V/R 
0  -  1/2  <v,x  -  u,y> 


C?> 


where  R  is  the  radius  of  curvature  of  the  panel  at  the  aid- 
thickness  plane.  These  values  for  the  strains  and 
curvatures  are  given  by  Sanders'  kinematic  relations 
considering  no  initial  shell  imperfections  C193. 

By  using  the  Kirchof f-Love  hypothesis  and  knowing  the 


Middle  surface  strains  and  curvatures,  the  strain-curvature 
relationship  for  a  lasinate  becoees 


By  substituting  the  strain  variation  through  the 
thickness.  Equation  (10),  into  the  stress-strain  relations. 
Equation  (6),  the  stress  in  the  layer  can  be  expressed 
in  teras  of  the  lasinate  Middle  surface  strains  and 
curvatures  C293  as 


share  the  elements  are  as  defined  in  Equation  (5> . 

The  resultant  forces  and  eoaents  acting  on  a  lasinate 


are  obtained  by  integrating  the  stresses  in  each  layer  of 
the  laainate.  For  as  a  force  per  unit  length  and  as  a 
aoaent  per  unit  length.  Me  have 


/ 1 

/: 


<TX  dt 


<rK  *  dz 


The  forces  and  moments  acting  on  a  flat  laminate  are 
shown  on  Figure  2.3. 
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For  an  N— layered  laminate. 


the  force  and 


t 


resultants  are  C263 


where  the  kth  and  k-1  layers  are  as  defined  in  Figure  2.2. 

If  one  substitutes  Equation  (11)  into  Equation  (13)  and 
takes  into  account  the  fact  that  the  stiffness  aatrix  is 
constant  within  the  lamina,  the  stiffness  aatrix  can  be 
pulled  outside  the  integration  over  each  area  and  then 
summed  with  the  resultants  for  each  layer.  The  force  and 
moment  resultant  equations  become 


Since  the  strains  and  curvatures  in  Equation  (14)  are 
middle  surface  values  and  not  functions  of  z,  they  can  be 
removed  from  under  the  summation  signs  so  that  Equation  (14) 
can  be  written  as 
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bending  stiffnesses. 
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The  basic  solution  technique  followed  is  to  use  the 
finite  eleeent  Method  in  order  to  find  the  displsceeent 
vector,  calculate  the  Middle  surface  strains  and  curvature 
using  Sanders’  equations,  calculate  the  transformed  reduced 
stiffness  Matrix  and  then  determine  the  resultant  laminate 
forces  and  moments.  All  the  applied  forces  are  known  so 
that  the  potential  energy,  as  used  in  the  finite  element 
method,  can  be  calculated  and  thus  a  check  for  static 
equilibrium  using  a  nonlinear  solution  algorithm  can  be 


carried  out 


STAGSC-1  Thtpry 


The  Structural  Analysis  -for  General  Shells  (STAGS) 
computer  program  is  an  energy  based  finite  element  program 
that  has  been  under  development  by  the  Lockheed  Palo  Alto 
Research  Lab  since  1967.  The  latest  version,  STAGSC-1,  is 
finite  element  based  and  uses  a  Newton— Raphson  nonlinear 
equation  solver.  STAGSC-1  is  a  general  purpose,  thin— shell, 
structural  analysis  program  that  is  designed  principally  for 
the  nonlinear  static  and  dynamic  analysis  of  thin  shells 
C27].  For  studies  of  the  collapse  of  shells  with  cutouts, 
the  critical  axial  load  is  determined  by  the  use  of  a 
nonlinear  analysis  C23. 

Sob el  C271  points  out  that  the  geometrically  nonlinear 
analysis  in  STAGSC-1  is  based  on  small  engineering  stress- 
strain  relations  («  <.1>,  moderate  rotations  (less  than  0.3 
radians)  and  incremental  solution  of  force-displacement 
equations  with  the  modified  Newton— Raphson  (MNR)  method. 
In  the  MNR  method  the  stiffness  matrix  is  factored  and  its 
inverse  is  used  instead  of  reforming  the  matrix  for  each 
load  step.  The  matrix  is  inverted  when  the  solution  is  slow 
to  converge  or  diverges.  The  stiffness  matrix  is 
periodically  updated  (refactored)  within  the  MNR  method. 
Due  to  the  nonlinearity  of  the  problem,  the  geometric 
stiffnesses  and  displacements  are  input  into  the  stiffness 
matrix.  At  each  load  increment,  a  new  set  of  displacement 
functions  is  developed.  Periodical ly,  the  stiffness  array 
is  updated  to  include  the  new  displacement  functions. 
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The  static,  nonlinear  p rob  1m  is  a  solution  of  the 
equations  of  equilibriue  with  a  Taylor  expansion  of  the 
total  force  vector  about  the  currently  deforeed  state. 
The  summation  of  forces  includes  all  applied,  restoring  and 
residual  force  vectors  present  in  the  static  equilibrium 
expressions  and  can  be  t«ritt»n  as  1273 


{Fix  )>  +  d<F> { (tx  -  (x  ))  ♦  H.O.T. 
n  ■■ :  n+i  n 

d(x)  ,  .  ,  n_ 

|tx>*tx  > 

«  F (x  )  +  d(F>  I  <4x>  ♦  H.O.T. 
d<K>  tx>«txn> 


where  H.O.T.  represents  higher  order  terms.  In  this 
notation.  Fix)  is  the  sum  of  the  applied,  restoring,  and 
residual  force  vectors.  The  value  tx>  defines  the  currently 

deformed  state.  Note  that  <6(x)>  I  means  the 

tx>*Oc  > 
n 

evaluation  of  the  vector  C6(x)>  at  the  values  of  Cx>. 
Therefore, 


— dCFix  )> 

- - -  l  Ax)  -  <F(x  >>  ♦  H.O.T.  C18) 

dfx  >  n 

n 

If  one  notes  that  the  derivative  of  the  force  vector  is  the 
negative  of  the  nonlinear  stiffness  matrix,  CK(xn>3, 
Equation  <18)  may  be  written  as 


CKix  ) 3  < 
n 

x 3  •  IF (x  >>  +  H.O.T.  -  0 
n 

(19) 

The 

STAGSC-l 

program  treats  the  left 

hand 

side  of 

Equation 

(19)  as 

a  nonlinear  operator  L 

acting 

on  the 
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incrwntal  displacwntt,  and  thus  Equation  (17)  wy  b* 
written  as 

.«>>  *  <L(x  >>  -  <R>  ♦  H.O.T.  -  0  (20) 

nt-i  n 

Mhara  (R)  is  the  vector  o-f  externally  applied  -forces  and  the 
higher  order  teres  are  generated  by  the  strain  energy  and 
residual  forces.  The  operator  L  is  the  first  derivative  of 
the  strain  energy  functional  C273.  The  MNR  algorithm  used 
in  STASSC— 1  can  be  written  as 

Cx  -  <x  >  «  CL*  (x  )l“1  (CR>  -  CL(x  )>)  (21 ) 

n+i  n  n  n 

where  CL* (x  )]  1  is  the  factored  stiffness  eatrix  at  soee 
n 

previous  step.  The  nonlinearities  included  in  the 
coeputation  of  L(xn)  are  purely  geometric.  Figure  2.4 
represents  the  MNR  algor it ha  used  in  STABSC-1  as  a 
conceptual  flowchart.  The  STA6SC-1  calculation  of  the 
vector  <L(x)>  is  analogous  to  the  Method  of  calculating  a 
pseudo-force  vector  to  account  for  nonlinearities  C27J. 

Since  the  sue  of  the  forces  is  equal  to  zero  (Equation 
<17) )f  the  total  potential  energy  nust  be  stationary.  The 
total  potential  energy  of  a  body  is  equal  to  its  internal 
strain  energy  ninus  the  product  of  the  external  forces  and 
deflections  and  can  be  written  as 

V  -  U  -  W  (22) 

where  N  is  the  work  done  and  U  (Equation  (23) )  is  the  strain 
energy  of  the  body. 


Thwfori 


C  a  €  *  o  €  +  o  y  >  d  Vol 

xx  v  v  xy  'icy 


Hi*  values  for  the  stresses  and  strains  are  calculated 


as  in  the  Classical  Laaination  Theory  section  of  this 
thesis. 

ST AG SC— 1  uses  a  Choi esky  triangular  decoeposition  with 
forward  and  backward  substitution  to  solve  the  equations. 
Finite  eleeent  techniques  are  used  to  foraulate  the 
displacement  vector  in  teres  of  the  eleaent  shape  functions 
and  nodal  degrees  of  freedoa.  The  displaceaent  vectors  are 


finally  calculated  using  a  aodified  Newton— Raphson  iteration 
scheae.  Once  this  displaceaent  vector  is  known,  stresses. 


strains  and  resultant  forces  can  be  deterained  for  the 
increeental  loadstep. 


the  finite  eleaent  analysis,  a  graphite  epoxy 


cooposite  panel  was  aodeled.  The  panel  was  12  inches  long 
with  a  chord  length  of  12  inches,  and  a  radius  of  curvature 
of  12  inches.  The  panel's  diaensions  and  aaterial 
properties  are  the  saae  as  those  used  by  Hebert  C221  and 
Janisse  C233  in  their  work.  Figure  2.S  shows  the  panel's 


aaterial  properties,  diaensions  and  sign  conventions.  The 
cutout  is  shown  on  Figure  2.S  with  diaensions  a  and  b.  The 
values  for  the  corresponding  diaensions  of  the  notches  and 
cutouts  studied  are  listed  in  Table  2.1.  The  panel  has  a 


'•Aw  .■  • '  -  • 
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Di Mnii on*  and  Material  Properties 
Hatarial:  Graph i t a— Epoxy 

Radius:  R  *  12" 

Langth:  L  ■  12" 

Nuabar  of  PI las:  8 

Orientation  of  Plies,  0:  (0, 45, -45, 90) s 

Thickness:  8  plies  at  .005"  -  .04" 

Elastic  Moduli:  Ef  -  20500  ksi,  E2  ■  1300  ksi 

Shear  aodulus:  6  ■  750  ksi 

Poisson’s  Ratio:  ■  .33 

x,  y,  z:  Structural  Coordinates 

u,  v,  m:  Displacements 

Ru,  Rv,  Rw:  Rotations 

Figure  2.5.  Panel  Notation  and  Material  Properties 


Table  2.1 


Di sansi one  of  Cutouts  Analyzed 

Size  of  Cutout 
(inches) 

a  b 

2  2 
1  2 
2  2 
2  4 

4  2 

4  4 

synnetric  ply  layup  E0,45,— 45,903g  consisting  of  8  plies  of 
0. 005  inches  each  for  a  total  panel  thickness  of  .04  inches. 

For  all  but  the  2"x2"  notch,  a  syaaetric  set  of 
boundary  conditions  were  used.  The  tap  of  the  panel  was 
claaped  with  only  the  u  di spl acseent  free.  The  vertical 
edges  of  the  panels  had  the  u  and  v  displaceeents  (x  and  y 
direction  displaceeents,  respectively)  and  rotation  about 
the  x  axis,  Ru,  free.  The  bottoe  of  the  panel  was  claeped 
with  no  free  degrees  of  fraedoe  (DOFs).  Refer  to  Figure  2.S 
for  displaceeent  and  rotation  vectors. 

The  2”x2”  notch  was  used  for  the  boundary  condition 
study  in  which  unsyeeetric  sets  of  boundary  conditions  were 
evaluated  in  addition  to  the  syemetric  boundary  conditions 
used  for  the  other  cutouts.  The  four  sets  of  unsyeeetric 
boundary  conditions  evaluated  are  listed  in  Table  2.2. 


Table  2.2 


Unsymctric  Boundary  Conditions  for  the  2"x2"  Cutout 


Boundary 

Condition 

1 

2 

3 

4 

Top  of 
panel 

u  free 

u  free 

u  free 

u  free 

Right  side 
of  panel 

u,  v,  Ru 

free 

Top  half: 
u, Ru  free) 
Bottom  half: 
u,  v,  Ru 

free 

u,  v,  Ru 

free 

u,  v,  Ru 

free 

Bottom  of 
panel 

clamped 

clamped 

clamped 

clamped 

Left  side 
of  panel 

Top  half 
u,Ru,free 
Bottom  half: 
u,  v,  Ru 

free 

u,  v,  Ru 

free 

Top  edge: 
u,  v  , Ru 
free 
Middle: 

Ru  free 
Bottom 
half: 

Top  1/4: 
u,  Ru 

free 

Bottom  3/4 
u,  v,  Ru 

free 

u,  v,  Ru 
f  ree 

ALL  OTHER  DOFs  WERE  FIXED 


eiBffgnt  Sritctign 


The  alaaant  salactad  for  this  analysis  is  a 
quadrilateral  plata  element  referred  to  as  the  QUAF  411  in 
the  STAGSC-1  documentation  C253.  This  is  the  same  element 
used  by  Janisse  C233.  The  QUAF  411  has  32  DOFs  I  3 
translational  rotations,  2  in  plane  rotations  and  2 
independent  normal  rotations  at  each  corner  node,  plus 
tangential  displacements  at  each  of  the  four  midside  nodes. 
Figure  2.6  shows  the  QUAF  411  with  its  DOFs.  The  QUAF  411 


V  allLSi 
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Figure  2.6.  Quadrilateral  411  Plate  Elenent  with 
Local  Coordinate  System 

element  was  developed  to  remove  the  displacement 
incompatibility  which  exists  when  flat  shell  elements  do  not 
lim  in  the  same  plane  C273.  If  the  edge  displacements  are 
to  be  compatible,  the  transverse  displacement  shape 
functions  must  be  of  the  same  order  as  the  in —plane  shape 
functions  which  are  cubic.  Element  411  has  two  rotations 
at  each  corner  node  which  permit  individual  rotation  of  each 
adjacent  side  and  therefore  allow  shear  strain  at  the 
corners.  The  shape  functions  for  the  411  element  are  cubic 


parallel  to  the  edge  and  quadratic  perpendicular  to  the  edge 


for  in-plan*  displacements.  The  bending  shape  function*  are 
cubic  in  both  directions  £273.  Since  the  panels  studied  in 
this  thesis  require  the  nonlinear  collapse  analysis  branch 
of  STAGSC-1  which  uses  large  displacement  theory  and 
moderate  rotations,  the  extra  DOF  at  the  corner  node  and  the 
translational  DOF  at  the  midside  node  sake  the  QUA*-  411 
element  the  best  available  element.  A  drawback  to  the  QUAF 
411  element  is  that  there  is  no  transition  element  in  the 
STA6SC— 1  library,  so  a  constant  grid  size  must  be  used  [271. 


& 


AND  RESULTS 


Finite  Element  Mesh 


A  mesh  refinement  study  mb  conducted  on  the  panel  with 
the  2"x2"  notch  to  determine  an  optimal  finite  element  mesh 


arrangement.  A  viable  tradeoff  between  improved  results  due 


to  an  increased  number  of  data  points  and  the  corresponding 
increase  in  computer  time  had  to  be  considered.  The  results 


from  this  convergence  study 


then  applied  to  the  design 


of  the  finite  element  mesh  for  other  panels. 

Initially,  a  finite  element  mesh  consisting  of  13  rows 


and  13  columns  of  one  inch  square  QUAF  411  elements  was 
tested.  This  mesh  is  referred  to  as  the  coarse  mesh.  In  a 
finite  element  method  for  solving  a  continuum  mechanics 
problem,  the  approximation  improves  as  more  elements  are 


used  C28I.  Therefore,  the  13x13  mesh  provides  the  roughest 
approximation  for  the  panel's  behavior.  The  coarse  mesh  was 
then  refined  by  cutting  the  elements  immediately  around  the 
cutout  into  one  half-inch  square  elements.  Due  to  the  fact 
that  there  is  no  transition  element  for  the  QUAF  411 
element,  some  elements  in  the  mesh  are  rectangular  when  the 
mesh  is  refined. 

The  rectangular  elements  have  aspect  ratios  that  are 


not  unity.  Aspect  ratio  will  be  defined  as  the  length  of 
the  element  in  the  x  (axial!  direction,  divided  by  the 
length  of  the  element  in  the  y  (circumferential >  direction. 
An  element  with  an  aspect  ratio  of  one  does  not  experience 


xv.v 


TABLE  3.1 

Mesh  Refinement  Results  for  a  2"  x  2"  Notch 

13x13  17x19  19x17  19x19  21x19  25x25 


(Active) 

Bandwidth 

CPU  Time  to 
Collapse 
Load  (sec) 


97 

129 

128 

145 

145 

192 

545 

1156 

925 

1097 

1582 

4852 

17.6 

208.1 

217.3 

207.3 

207.6 

191 

12.2 

8.2 

12.1 

7.9 

8.0 

- 

28 

26 

21 

20 

26 

22 

Collapse 

Load 

(lbs/in) 

Difference 
Froe  25x25 
6rid  (X) 

Number  of 
Load  Steps 


the  biases  introduced  when  using  eleeents  with  different 
Magnitudes  of  length  in  the  x  and  y  directions.  Cook  C283 
states  that  a  desirable,  but  not  Mandatory,  requirement  for 
convergence  is  that  an  eleaent  should  have  no  preferred 
directions.  In  other  words,  the  element  should  be 
geometrically  invariant.  Various  mesh  refinement 
arrangements  were  tested  to  determine  how  elements  with 
aspect  ratios  other  than  unity  affected  the  behavior  of  the 
panel . 

The  results  of  the  mesh  refinement  study  can  be  found 
in  Table  3.1.  In  this  table,  CPU  (Central  Processing  Units) 
time  is  the  amount  of  computer  time  used  to  reach  collapse 
load.  The  bandwidth  is  the  full  bandwidth  used  for  matrix 
storage  in  the  finite  element  analysis.  The  number  of  load 


steps  is  the  nuabar  of  tines  a  load  increeent  Mas  applied. 
As  collapse  load  is  approached,  the  computer  program  cuts 
the  load  increment  to  achieve  convergence.  The  results,  as 
shomn  in  this  table,  -for  the  25x25  grid  (consisting  entirely 
of  one  half— inch  square  elements)  are  taken  from  Janisse 
C233.  The  coarse  grid  had  a  collapse  load  that  mas  12. 2% 
higher  than  the  25x25  grid  but  took  88%  less  computer  time. 
Janisse  C233  also  found  that  the  collapse  load  for  a  coarse 
grid  Mas  18%  higher  than  experimental  results,  referred  to 
subsequently,  (experimental  collapse  load  Mas  2121  lbs) 
Mhereas  the  25x25  mesh  had  a  collapse  load  that  Mas  only 
7.5%  higher  than  the  experimental  collapse  load.  The  mesh 
refinement  study  is  a  method  for  determining  the  best  May  to 
more  closely  approximate  experimental  results,  Mithout  using 
an  unreasonable  amount  of  computer  time. 

Figure  3.1  shoes  the  refined  mesh  in  the  vicinity  of 
the  cutout  for  each  of  the  different  mesh  arrangements 
studied.  By  comparing  the  results  for  the  17x19  and  19x17 
grid,  one  can  observe  that  the  17x19  grid  more  closely 
approximates  experimental  collapse  load  but  uses  20%  more 
computer  time  than  the  coarse  mesh.  The  19x17  mesh  does  not 
shorn  an  improvement  in  collapse  load,  indicating  that 
elements  Mith  an  aspect  ratio  of  .5  hinder  the  advantages 
that  are  expected  Mhen  refining  a  mesh.  The  19x19  mesh  had 
half-inch  square  elements  a  uniform  distance  of  2  inches 
from  the  opening.  This  refinement  Mas  found  to  more  closely 
approximate  the  collapse  load  Mhile  using  23%  less  computer 
time  than  the  25x25  grid.  A  further  grid  refinement  to  a 
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21x19  Msh  did  not  iaprovt  the  col  laps*  load  approximation 
over  values  obtained  for  the  19x19  mesh,  but  used  aleost  900 
seconds  aore  computer  time  than  the  19x19  mesh.  From  a 
collapse  load  approximation  point  of  view,  the  19x19  mesh  is 
the  best  refined  mesh  to  choose. 


To 

f urther 

analyze 

the  meshe 

s,  the  values 

for 

the 

moments. 

"v 

and 

Mx,  «nd 

force  ,N 

y,  were  compared 

in 

the 

vicinity 

of 

the 

notch. 

(Refer 

to  Figure  2.3 

for 

the 

individual  forces  and  moments  acting  on  a  laminate).  These 
force  parameters  mere  chosen  for  analysis  because  they  are 
important  components  in  a  collapse  analysis.  Each 
functional  value  was  compared  along  rays  at  different  angles 
extending  from  the  upper  left  corner  of  the  notch.  Figure 
3.2  is  a  comparision  of  versus  distance  from  the  corner 
of  the  notch  for  all  five  meshes  along  an  outward  ray  at 
135°  from  the  upper  left  corner  of  the  notch.  The  13x13 
mesh  provides  the  least  accurate  approximation  for  behavior 
because  it  has  fewer  elements.  As  the  mesh  is  refined,  the 
values  of  for  each  of  the  cases  tested  begin  to  cluster 
together.  The  19x17  mesh  does  not  follow  the  convergence 
pattern  of  the  other  three  refined  meshes.  Over  1.5  inches 
from  the  opening  where  the  rectangular  elements  start,  the 
values  for  My  diverge.  The  elements  with  a  .5  aspect  ratio 
are  only  an  inch  away  from  the  opening  thus  not  allowing  an 
improvement  in  results  as  would  be  expected  with  a  refined 
mesh.  The  .5  aspect  ratio  elements  are  longer  in  the 
circumferential  direction  and  therefore  the  moment  about  the 
x  axis,  M ,  appears  to  be  adversely  affected  by  the 
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Figure  3.2.  Moment  about  X  Axis  versus  Distance  from 
Corner  along  a  Ray  at  135° 
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shorter  dimension  in  the  axial  direction.  An  eleeent  with 

an  aspect  ratio  of  1  does  not  seee  to  experience  the 

difficulties  associated  with  an  oblong  eleeent.  The  results 

for  the  19x19  and  21x19  aesh  do  not  differ  significantly. 

To  exaeine  the  eeshes  further,  values  for  along  a 

line  at  45°  froe  the  corner  are  coapared  in  Figure  3.3. 

Even  though  the  aagnitude  of  the  eoeent  is  saall,  the  values 

for  each  of  the  aeshes  are  consistent  froa  a  finite  eleeent 

aodel  point  of  view.  As  the  aesh  is  refined,  the  values  for 

Mk  begin  to  cluster  together  except  for  the  coarse  aesh. 

The  17x19  aesh  diverges  froa  the  others  over  l.S  inches  froa 

the  corner  where  eleaents  with  an  aspect  ratio  of  2.0  start 

appearing  in  the  aesh.  The  divergence  for  the  17x19  aesh  is 

noted  in  Figure  3.3  whereas  it  did  not  happen  in  Figure  3.2. 

This  is  due  to  the  fact  that  Figure  3.3  is  a  plot  of  aoaent 

about  the  y  axis,  H^,  so  eleaents  with  the  longer  diaension 

in  the  vertical  direction  affect  the  behavior  of  the  panel. 

This  leads  to  the  conclusion  that  the  aesh  should  be  refined 

equidistant  around  the  cutout,  as  can  be  seen  froa  the  19x19 

aesh  arrangeaent  in  Figure  3.1. 

A  further  refineaent,  the  21x19  grid,  does  not  iaprove 

the  results  for  M  or  M  over  the  19x19  aesh.  The  21x19 

x  y 

also  does  not  adversely  affect  the  values  for  M  or  H  . 

x  y 

This  is  because  the  oblong  eleaents  are  far  enough  away  froa 
the  cutout  so  as  not  to  cause  problems  with  the  results. 
Since  the  21x19  aesh  does  not  iaprove  results  but  uses  no re 
coaputer  tiae,  the  19x19  aesh  appears  to  be  the  optiaua 
finite  eleeent  aesh  for  the  2"x2"  notch. 
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Figure  3.3.  Moment  about  Y  Axie  versus  Distance  from 

o 

Corner  along  a  Ray  at  45 
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One  additional  paraactar,  fore*  in  tha  y  direction,  N  , 
mi  examined  (Figure  3.4)  along  a  ray  at  90°.  The  eemh 
characteristics  follow  the  ref insesnt  in  a  sieilar  fashion 
to  the  aooent  function  coeparisons.  All  of  the  refined 
eeshes  are  good  approximations  over  the  coarse  mesh  for  the 
force  function. 

From  examining  the  mesh  parameters,  the  13x13  mesh  is 
the  least  accurate  due  to  the  fewer  number  of  elements  used 
in  the  numerical  approximation.  With  the  17x19  and  the 
19x17  meshes,  difficulties  arise  when  considering  the 
moments  about  axes  in  the  direction  of  the  shorter  dimension 
of  the  oblong  element.  The  19x19  mesh  has  elements  with  an 
aspect  ratio  of  1  over  a  uniform  distance  from  the  opening. 
This  mesh  experienced  no  irregularities  when  approximating 
moments.  A  more  refined  mesh  (21x19)  gave  no  significant 
improvement  for  values  in  the  vicinity  of  the  notch.  Oblong 
elements  greater  than  2  inches  from  the  cutout  had  no 
adverse  effects  on  the  buckling  behavior  of  the  panel. 

The  results  from  the  mesh  refinement  study  were  next 
applied  to  the  finite  element  mesh  for  the  other  panels 
tested.  As  stated  earlier,  it  is  best  to  refine  the  mesh 
evenly  around  the  cutout.  If  this  is  not  possible  due  to 
the  size  or  positioning  of  the  cutout,  elements  with  an 
aspect  ratio  other  than  unity  should  be  at  least  four  square 
elements  from  the  cutout.  Elements  with  aspect  ratios  of  .5 
adversely  affect  and  collapse  load,  and  should  be  at 
least  2  inches  from  the  cutout.  Elements  with  aspect  ratios 
of  2  affect  and  should  be  at  least  i.S  inches  from  the 
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Figure  3.4. 


Axial  Force,  Nyt  versus  Distance  f roe 

o 

Corner  along  a  Ray  at  90 


cutout 


The  Mih  arrangi 


its  chosen  for  each  of  the 


notchas  will  ba  discussed  latar  in  thia  thaaia. 

Bgyntf^ry  C.g04iU0h  Study 

An  axpariaantal  study  has  previously  baan  carried  exit 
by  the  Fatigue,  Fracture  and  Reliability  Group  of  the  Air 
Force  Wright  Aeronautical  Laboratories  at  Wright 'Patterson 
AFB  considering  a  2"x2"  notch  within  the  sane  type  of 
grapite-epoxy  panel  described  previously.  An  atteept  to 
siaulate  a  sinple  support  boundary  condition  along  the 
vertical  edges,  and  claaped  at  the  tap  and  botton  was 
carried  out.  Results  f  roe  this  study,  referred  to  in 
subsequent  paragraphs,  indicated  that  the  experiaental 
boundaries  were  not  acting  as  originally  assuaed,  yet 
collapse  load  was  within  approxinately  10%  of  the 
experiaental  value.  Thus,  the  author  conducted  a  finite 
eleaent  study  of  varieties  of  boundaries  to  try  to 
approxiaate  the  actual  boundaries  present  within  the 
experiaent.  It  is,  to  the  best  of  the  author's  knowledge, 
the  first  atteept  at  nunerically  siaulating  an  experiaental 
*et  of  boundaries  for  a  composite  panel  with  a  given  size 
notch. 

Several  different  boundary  condition  cases  were  studied 
for  the  2"x2“  notch  with  the  refined  aesh  in  an  attempt  to 
•ore  realistically  model  the  experieental  results.  Four 
different  sets  of  unsyaaetric  conditions  were  modeled  (Table 
2.2)  along  with  the  symmetric  boundary  conditions  described 


Figure  3.5.  Test  Apparatus  Mith  Panel 


in  Section  II  of  this  thesis. 

The  experiaental  text  fixture  (Figure  3.5)  Mas  siailar 
to  that  used  by  Hebert  C22]  and  Janisse  C233.  The  panel  Mas 
fabricated  from  NARMCO  T300/5208  graphite-epoxy.  A  two  inch 
square  hole  Mas  cut  out  of  the  center  of  the  panel.  The 
corners  of  the  cutout  Mere  aachined  to  be  slightly  rounded 
<1/16  inch  radius)  rather  than  a  perfect  square  cutout  so 
that  there  are  no  singularities  at  the  corners.  For  a  eore 
detailed  description  of  the  test  apparatus,  see  Reference 
C223.  The  experimental  apparatus  Mas  modified  for  this 
testing  by  placing  Linear  Variable  Differential  Transformers 
(LVDTs)  at  various  locations  on  the  panel  to  measure  out— of- 
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Figure  3.6.  LVDT  Location  and  Initial  Radial  Diaplacaaant 
on  Experimental  Panel 
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plan*  (radial)  displacements.  Figure  3.6  shows  the 
locations  of  the  LVDTs  on  the  panel  as  viewed  from  the  front 
of  the  panel  (panel  radius  into  the  page) . 

Various  analytic  boundary  conditions  were  examined  in 
addition  to  the  symmetric  boundary  conditions,  because 
previous  studies  did  not  model  experimental  results 
precisely.  During  the  experimental  testing,  depressions  in 
the  teflon  tape  that  held  the  panel  into  the  side  restraints 
were  noted.  The  teflon  tape  on  the  left  side  of  the  panel 
had  the  largest  depressions.  This  indicated  that  the 
experimental  boundary  conditions  were  changing  from  the 
initial  conditions.  An  attempt  at  modelling  this  phenomena, 
was  carried  out  by  fixing  the  v  displacement  first  on  the 
upper  half  of  the  left  side  of  the  panel,  and  then  on  the 
upper  half  of  the  right  side.  Next,  the  v  displacement  was 
held  on  the  upper  quarter  of  the  left  side  of  the  panel. 
Finally,  both  the  u  and  v  displacements  were  fixed  between  1 
inch  to  6  inches  down  the  left  edge  of  the  panel.  These 
results,  along  with  the  results  for  the  symmetric  boundary, 
were  compared  to  the  experimental  findings.  Radial 
displacements  at  several  different  LVDT  locations  were 
compared  as  a  function  of  load. 

The  experimental  panel  was  also  examined  in  order  to 
see  if  there  was  any  initial  eccentricity  or  displacements 
introduced  into  the  panel  by  the  test  apparatus  prior  to  the 
loading.  The  initial  radial  displacements  of  the  panel  were 
measured  by  us*  of  a  set  of  dial  vernier  calipers.  This 
device  was  also  used  by  Hebert  [22]  to  measure  surface 
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imperfections.  Using  0.0  as  a  datum  surface*  positive 
displacement  is  measured  radially  outward.  The  top  and 
bottom  of  the  panel  were  designed  to  be  at  a  displacement  of 
0.0.  Theoretically*  if  the  panel  is  manufactured  and 
mounted  perfectly*  the  entire  panel  would  have  an  initial 
radial  displacement  of  0.0.  Figure  3.6  shows  the  initial 
radial  displacements  at  each  of  the  LVDT  locations. 

At  the  point  corresponding  to  LVDTs  2  and  17  there  is 
an  outward  bulge*  whereas  in  the  vicinity  of  the  cutout* 
there  is  a  negative  displacement.  These  eccentricities  can 
mainly  be  attributed  to  initial  inaccuracies  introduced  by 
the  test  facility  when  mounting  the  panel  in  the  test 
apparatus.  The  effect  that  these  initial  eccentricities 
have  on  the  actual  panel  as  compared  to  analytic  results  can 
be  seen  in  Figures  3.7  through  3.11.  In  these  figures* 
boundary  conditions  1,  2*  and  3  are  as  listed  in  Table  2.2. 
Boundary  condition  4  in  Table  2.2  was  considered  but  not 
plotted  because  the  results  did  not  differ  significantly 
from  the  results  for  boundary  condition  1. 

Analytical  and  experimental  results  for  the  radial 
displacement  at  LVDT  location  2  are  compared  in  Figure  3.7. 
One  can  observe  that  the  variety  of  boundary  conditions 
tested  fall  significantly  away  from  experimental  data.  The 
author  notes  that  the  boundary  condition  in  itself  does  not 
solely  affect  the  displacement  of  the  panel.  The  effect 
that  the  overall  initial  negative  displacement  has  on  the 
experimental  results  at  LVDT  location  2  can  be  seen  in 
Figure  3.7.  The  experimental  data  shows  that  the  panel  at 
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point  2  is  initially  displaced  outward*  but  as  tha  load  is 
applied*  tha  point  starts  to  bow  out  and  than  appaars  to 
raversa  its  aotion  and  aaintains  an  alaost  linear  relation 
as  a  -function  of  load.  A  possible  course  -for  this  phenomena 
aay  be  explained  by  the  fact  that  the  negatively  displaced 
area  around  the  cutout  acts  as  a  restraining  force  analagous 
to  a  support.  The  area  absorbs  moments  and  acts  siailar  to 
a  depressed  clamped  edge.  This  appears  to  cause  a  constant 
displacement  function  of  experimental  data  for  point  2 
because  the  cutout  area  is  holding  the  panel  and  preventing 
it  from  moving  in  the  direction  of  its  initial  displacement. 
The  analytic  results  for  all  of  the  boundary  conditions 
tested  are  very  similar  and  do  not  approximate  the 
experimental  results  at  point  2  as  well  as  one  would  like. 

By  comparing  analytic  and  experimental  results  for 
radial  displacement  at  LVDT  location  3*  Figure  3.8*  one  can 
see  that  the  results  for  boundary  condition  3  better 
approximate  with  experimental  data.  At  point  3*  the  results 
for  boundary  condition  3  show  almost  no  radial  displacement. 
This  implies  that  the  effect  of  boundary  condition  3  is  to 
physically  create  a  counterflexure  point  at  the  given  LVDT 
location  and  therefore  will  theoretically  have  no  moment  in 
the  y  direction.  The  effect  of  the  eccentricities  can  also 
be  observed  in  Figure  3.8.  The  cutout  area  is  more 
depressed  than  the  panel  around  point  3  and  therefore 
restrains  the  point.  tie  would  expect  point  3  to  move 
inward*  but  the  panel  in  the  vicinity  of  the  cutout 
restrains  the  movement  and  as  the  applied  load  reaches  500 
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Figure  3.7.  Radial  Displacement  Plot,  for  Boundary 
Conditions  Studied,  at  LVDT  Location  2 


40 


COMPRESSIVE  LORD  (LBS) 

0.0  500.0  1000.0  1500.0  2000.0  2500.0  3000. 


Figure  3.8.  Radial  Displacement  Plot,  for  Boundary 


Conditions  Studied,  at  LVDT  Location  3 
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lbs,  th*  displacement  reverses  since  the  area  around  the 
cutout  is  acting  analogous  to  like  a  support.  Zt  therefore 
appears  that  point  3  is  experi mentally  undergoing 
displacement  as  a  function  of  surface  imperfections  and  a 
change  in  boundary  condition  similar  to  analytic  bounary 
condition  3  which  is  a  restrained  edge. 

If  one  examines  the  panel  at  LVDT  4  in  Figure  3.9,  it 
is  observed  that  the  large  initial  displacements  at  the 
upper  right  corner  of  the  cutout  (corresponding  to  LVDT  6> 
acts  analogous  to  a  folded  plate  and  causes  point  4  to 
displace  outward.  There  is  a  tendency  for  the  area  of  the 
cutout  to  restrain  areas  of  the  panel.  Although  there  is  no 
applied  force  at  the  cutout,  the  restraint  generated  is  due 
to  the  radial  location  of  the  cutout.  The  area  of  the 
cutout  is  initially  displaced  inward  more  than  LVDT  location 
4.  This  causes  more  moment  restraint  to  act  on  point  4  with 
increasing  load  as  the  area  around  the  cutout  acts  similar 
to  a  folded  plate.  This  phenomena  can  only  be  modeled 
analytically  if  the  surface  eccentricity  is  precisely 

incorporated  into  the  finite  element  model.  These 

displacement  characteristics  appear  to  be  happening 
experimentally.  At  an  applied  load  of  about  1500  lbs,  the 
experimental  data  shows  that  point  4  starts  shedding  or 
redistributing  its  load  and  begins  to  move  radially  inward. 
Boundary  conditions  1,  2,  and  the  symmetric  conditions 

approximate  the  behavior  well  up  to  about  1500  lbs. 
Boundary  condition  3  bounds  the  experimental  data  on  the 
other  side.  This  seems  to  indicate  the  boundary  conditions 
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Figure  3.9.  Radial  Displacement  Plot,  for  Boundary 
Conditions  Studied,  at  LVDT  Location  4 
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tasted  in  this  thesis  for  point  4  bound  the  experimental 
displacement  data,  and  the  actual  conditions  to  be  applied 
analytically  are  somewhere  in  between  these  cases.  This  is 
due  to  the  fact  that  boundary  condition  3  was  applied  from 
the  onset  of  the  load  application  process.  The  experimental 
results  indicate  that  there  is  some  sort  of  restraint 
affecting  point  4.  Boundary  condition  3  best  models  this  at 
higher  loads  when  the  change  begins  to  occur  experimentally. 

Figure  3.10  compares  the  analytic  and  experimental 
results  for  radial  displacement  at  LVDT  location  11. 
Boundary  condition  3  generally  shows  the  same  trend  of 
displacement  as  experimental  results.  When  the  effect  of 
the  boundary  condition  is  considered  along  with  initial 
surface  imperfections,  the  translation  between  the  two 
curves  can  be  attributed  to  the  effect  of  the  initial 
eccentricities. 

Boundary  condition  3  does  not  affect  the  results  for 
displacement  at  point  17  (Figure  3.11)  as  it  did  for  the 
other  points  considered.  This  is  largely  due  to  the  fact 
that  point  17  is  a  greater  distance  from  the  changing 
boundary  than  the  other  points.  All  the  analytic  boundaries 
approximate  the  radial  displacement  at  point  17  to  the  same 
degree  because  the  point  is  so  far  removed  from  the  boundary 
that  the  action  of  the  restraint  is  not  felt. 

In  addition  to  the  initial  eccentricities  causing  the 
test  panel  to  behave  differently  from  analytic  results,  it 
is  noted  that  each  point  on  the  test  panel  does  not  start  at 
a  datum  of  0.0,  but  rather  at  an  amount  of  the  initial 
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Figure  3.11.  Radial  Displacement  Plot,  for  Boundary 


Conditions  Studied,  at  LVDT  Location  17 


displacement  from  as  the  datum.  This  means  that  each  of  the 
experimental  curves  on  Figures  3.7  through  3.11  should  be 
shifted  by  the  amount  of  the  initial  displacement  at  each  of 
the  respective  LVDTs.  This  mould  mean  that  the  analytic 
results  for  boundary  condition  3  better  approximate 
experimental  displacement  results. 

Boundary  condition  3  attempts  to  model  some  friction 
appearing  in  the  v  direction  depicted  by  small  depressions 
registering  on  the  teflon  tape  in  the  test  apparatus. 
Analytically,  a  portion  of  the  left  side  of  the  panel  was 
held  rigidly  with  rotation  about  u  permitted.  The  right 
side  was  free  to  move  in  the  analytic  simulation.  The 
results  for  LVDTs  3  and  4  (using  boundary  condition  3) 
compared  more  favorably  because  the  left  edge  boundary 
condition  more  accurately  models  what  was  happening 
experimentally.  There  is  more  displacement  allowed  further 
from  the  simulated  rigid  body  than  near  it  (LVDT  2).  A 
clamped  boundary  acts  as  a  fold  which  produces  bending 
moment  as  load  increases,  thus  causing  smaller  displacements 
near  the  edge.  As  the  applied  load  increases 
experimentally,  the  radial  displacements  become  large  and 
the  boundary  conditions  change  as  the  panel  buckles  outward. 
Portions  of  the  panel  surface  come  in  contact  with  the  side 
restraints  causing  friction.  The  results  at  LVDTs  11  and  17 
are  not  affected  as  dramatically  because  their  location  is 
farther  away  from  the  changing  boundary.  Boundary  condition 
3  is  not  a  true  approximation  of  the  experimental  boundary 
but  rather  appears  to  be  an  upper  bound.  The  true  boundary 


is  s  function  of  load  and  di spl acsaant .  To  battar  aodal 
experimental  rasults  analytical iy,  this  claapad  boundary 
condition  should  not  ba  appliad  at  tha  baginning  of  loading 
but  rathar  whan  tha  displacaaants  start  getting  large  and 
tha  depressions  fora  in  tha  teflon  tape. 

One  difference  that  occurred  Mith  boundary  condition  3 
is  that  tha  collapse  load  Mas  16X  greater  than  for  the  other 
conditions  tested  (Figure  3.12).  Boundary  condition  3  has  a 
collapse  load  of  greater  than  3000  lbs  Mhereas  the  other 
boundaries  tested  had  collapse  loads  of  approximately  2S00 
lbs.  The  higher  collapse  load  is  due  to  boundary  condition 
3  producing  a  greater  stiffness  (slope  of  load  versus  end 
shortening  curve)  in  the  panel  by  restraining  the  u 
displaceaent.  This  stiffness*  Mhich  is  tMice  that  for  the 
other  boundary  conditions*  leads  to  a  17X  higher  collapse 
load.  The  experimental  panel  has  initial  eccentricities 
Mhich  result  in  a  collapse  load  of  only  2200  lbs.  The 
geometrically  perfect  panels  used  in  the  analytical  Mork 
therefore  have  consistently  higher  collapse  loads  than  the 
experimental  panel.  Another  reason  that  the  collapse  load 
for  boundary  condition  3  is  higher  than  experimental  results 
is  that  analytically  the  u  restraint  Mas  applied  at  the 
beginning  of  the  load  application  process.  Experimentally* 
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representation.  In  fact,  the  boundary  condition  which  beet 
represented  the  experimental  displacement  yielded  a  collapse 
load  approximately  30%  greater  than  the  test  value.  This 
appears  to  indicate  the  fact  that  collapse  is,  considering 
this  particular  geometric  case,  a  function  of  axial  load 
which  is  comparatively  independent  of  radial  displacement. 
Thus,  one  appreciates  the  statement  that  buckling  is  an 
averaging  phenomena  in  which  most  secondary  effects  are 
unimportant.  Of  course,  it  is  required  to  keep  in  mind  that 
the  problem  just  described  relates  to  a  small  cutout  within 
a  cylindrical  panel.  The  relative  movement  within  the  panel 
is  primarily  linear.  A  subsequent  section  (Notch  and  Cutout 
Study)  of  this  thesis,  considers  several  larger  cutouts 
which  yield  truely  nonlinear  phenomena. 

Discussion  up  to  this  point  has  been  involved  with 
comparing  analytic  and  experimental  results.  Now,  this 
thesis  will  deal  with  a  strict  analytic  comparison.  Moments 
and  displacements  will  be  examined  in  order  to  gain  a  better 
understanding  of  the  behavior  of  a  cylindrical  panel  with 
the  various  boundary  conditions  previously  compared.  The 
radial  displacement  was  plotted  as  a  function  of  distance 
from  the  left  edge  of  the  panel  along  two  lines,  3  inches 
and  9  inches  from  the  top  of  the  panel.  The  displacement 
along  each  of  these  lines  was  plotted  for  a  low  load  level 
and  also  for  a  collapse  load  for  each  of  the  boundary 
conditions  tested  (Figures  3.13  through  3.16).  From  Figure 
3.13  it  can  be  seen  that  the  top  of  the  panel  deflects  more 
near  the  left  edge  of  the  top  of  the  panel  than  it  does 
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Not*:  Displacements  taken  along  lines  3”  and  9"  from 
top  edge  of  panel. 


Figure  3.13.  Radial  Displaceeent  as  a  Function  of  Panel 
Location  for  Symmetric  Boundary  Conditions 
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•long  the  line  9  inches  down.  There  is  a  trough  forcing 
radially  inward  from  the  upper  left  edge  of  the  cutout  to 
below  the  lower  right  edge  of  the  cutout.  Figures  3.14  and 
3. IS  show  that  when  v  is  restrained  on  one  edge,  the  radial 
displaceeent  is  not  affected  significantly  except  to  lessen 
the  radial  displaceeent  slightly  near  the  restrained  edge. 
By  restraining  u  displacement  (boundary  condition  3)  the 
radial  displacement  is  less  near  the  left  edge  as  cam  be 
seen  in  Figure  3.16.  The  trough  of  radial  displacaMaent  has 
shifted  dramatically  when  compared  to  the  other  three 
boundaries  studied.  By  restraining  the  u  displacement,  the 
trough  starts  above  the  upper  right  edge  of  the  cutout  and 
extends  diagonally  to  below  the  lower  left  edge  of  the 
cutout. 

The  values  for  M  and  H  ,  for  each  of  the  boundary 

1C  jr 

conditions  were  compared  along  lines  2.5  inches  and  9.5 
inches  from  the  top  of  the  panel.  The  moment  about  the  y 
axis  for  symmetric  boundary  conditions  (Figure  3.17), 
boundary  condition  1  (Figure  3.18),  and  boundary  condition  2 
(Figure  3.19)  are  very  similar,  and  have  no  significant 
differences.  The  maximum  and  minimum  moments  correspond  to 
the  points  of  greatest  outward  or  inward  radial 
displacement.  Restraining  the  v  displacement  on  the  upper 
half  of  the  panel  does  not  affect  significantly. 
Examining  for  boundary  condition  3  (Figure  3.20)  shows 
that  the  moment  is  smaller  for  the  2946  lb  load  near  the 
semi-clamped  left  edge.  The  moments  near  the  center  of  the 
panel  are  nearly  the  same  in  magnitude  even  though  boundary 
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Not*:  Displacement*  taken  along  lines  3"  and  9"  fr 
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Figure  3.14.  Radial  Displacement  as  a  Function  of  Panel 
Location  for  Boundary  Condition  1 
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Figure  3.17.  Moment  about  Y  Axis,  „  versus  Distanc 
from  Left  Edge  of  Panel  for  Symmetric 
Boundary  Conditions 
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Figure  3.18.  Moment  about  Y  Axis,  M^,  versus  Distance 
from  Left  Edge  for  Boundary  Condition  I 
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Figure  3.19.  Moment  about  Y  Axle,  M^,  versus  Distance 
from  Left  Edge  for  Boundary  Condition  2 
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Flgur*  3.25.  Mo**nt  about  Y  Axis,  «x,  versus  Distance 
fro*  Left  Edge  for  Boundary  Condition  3 


condition  3  had  a  higher  collapse  load  than  the  other  cases. 
This  is  due  to  the  restraint  on  the  left  edge  Mhich 
correspondingly  reduces  near  the  panel’s  center.  The 
values  for  near  the  right  edge  of  the  panel  for  boundary 
condition  3  are  greater  than  the  corresponding  eoeents  for 
the  other  boundary  conditions.  This  is  caused  by  the  right 
edge  being  freer  to  eove  as  coepared  to  the  left  edge.  The 
■aments  at  collapse  have  the  same  trend  for  the  symmetric 
boundary  and  boundaries  1  and  2  because  the  moments  are  not 
large  enough  to  be  affected  by  the  restraint  in  the  v 
direction.  Boundary  condition  3  produces  moments  that  are 
negative  at  the  left  edge. 

The  moment  about  the  x  axis  was  also  compared  along  the 
upper  line  (Figures  3.21  through  3.24).  Once  again  there  is 
no  significant  effect  introduced  by  restraining  the  v 
displacement  on  the  upper  half  of  the  panel.  For  boundary 
condition  3,  Figure  3.24,  the  values  for  are  much  less 
than  for  the  other  boundaries.  The  negative  values  for  My 
are  also  much  greater  between  7  and  9  inches  than  for  the 
previous  boundary  conditions.  This  occurs  because  the 
restrained  displacements  act  similar  to  a  clamp  such  that 
displacements  farther  away  oscillate  more  freely,  yielding 
larger  negative  moments  about  the  x  axis.  The  large  peak  in 
My  in  Figures  3.21  through  3.23  at  2.5  inches  from  the  left 
edge  corresponds  to  LVDT  location  2.  Boundary  condition  3 
more  accurately  models  experimental  results,  yet  there 
should  be  physically  less  moment  build  up  near  the  left  edge 
in  the  experimental  panel  than  the  amount  computed 
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Figure  3.21.  Moment  about  X  Axis,  M^t  versus  Distance 
from  Left  Edge  of  Panel  for  Symmetric 
Boundary  Condition 
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Figure  3.24.  Moment  about  X  Axis,  My,  versus  Distance 
from  Left  Edge  for  Boundary  Condition  3 
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analytically  becauu  the  axpariaantal  boundary  condition 
changes  with  applied  load.  In  each  of  these  figuresv  the 
values  for  eoeent  in  the  y  direction  go  to  zero  at  the  edges 
of  the  panel. 

The  values  for  H  and  M  were  also  coepared  along  a 

x  y 

line  9.5  inches  froa  the  top  of  the  panel.  This  line  is  far 
enough  away  from  the  changing  boundaries  so  that  one  would 
expect  little  change  between  the  values.  for  the 
symmetric  boundary  condition  (Figure  3.25)  does  not  differ 
significantly  froa  the  values  for  M^  with  boundary  condition 
1  (Figure  3.26)  or  2  (Figure  3.27).  These  two  unsyaaetric 
boundary  conditions  show  a  slight  increase  in  aoaent  at 
collapse  load  10  inches  froa  the  left  edge  of  the  panel.  A 
slight  influence  froa  the  restained  u  displaceaent  is 
experienced  by  M^  values  along  the  lower  line  for  boundary 
condition  3  (Figure  3.28).  The  aoaent  has  peaks  at 
approxiaately  the  same  locations  as  with  syaaetric  boundary 
condi tons,  but  the  magnitudes  are  slightly  less. 

Moment  about  the  x  axis,  M^,  was  also  exaai ned  at  the 
three  load  levels  along  the  lower  line.  Figure  3.29  shows 
the  results  for  the  symmetric  boundary  conditions. 
Comparing  boundary  condition  1  (Figure  3.30)  to  these  values 
shows  an  increase  in  M^  near  the  right  edge  of  the  panel. 
This  is  due  to  the  fact  that  the  entire  right  edge  of  the 
panel  is  freer  to  move  in  the  y  direction  as  compared  to  the 
left  edge  of  the  panel.  Figure  3.31  shows  this  same  result 
for  boundary  condition  2  except  the  phenomena  is  observed 
along  the  left  edge  of  the  panel.  The  magnitudes  for  M  are 


MX  (IN-LBS/IN) 


LP 


Not®:  Moments  taken  along  a  line  9.5"  f roe 

top  edge  of  panel. 


Figure  3.25.  Moment  about  Y  Axis,  M^#  versus  Distance 
from  Left  Edge  of  Panel  for  Symmetric 
Boundary  Condition 
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Figure  3.26.  Moment  about  Y  Axis,  M^,  versus  Distance 
from  Left  Edge  for  Boundary  Condition  1 
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Note:  Moments  taken  along  a  line  9.5“  from 
top  edge  of  panel. 

Figure  3.27.  Moment  about  Y  Axis,  M^,  versus  Distance 
from  Left  Edge  for  Boundary  Condition  2 


I-OJ 


MY  (IN-LBS/IN) 


Note:  Moments  taken  along  a  line  9.5"  from 

top  edge  of  panel. 


Figure  3.29.  Moment  about  X  Axis,  versus  Distance 
fro*  Left  Edge  of  Panel  for  Symmetric 
Boundary  Condition 
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grutar  for  boundary  condition  3  (Figura  3.32)  ax  cap  t  along 
tha  right  adga  of  tha  panal  where  thay  ara  slightly  lass. 

The  large  peaks  in  the  aoaant  about  tha  y  axis  along 
both  lines  near  the  aiddla  of  tha  panal  correspond  to  tha 
location  of  tha  cutout.  Along  tha  upper  line,  tha  aoaents 
about  tha  y  axis  are  negative  near  tha  canter  of  tha  panel, 
whereas  the  aoaents  are  positive  at  the  center  of  the  panel 
along  the  bottoa  line.  The  aoaents  are  syaaetric  about  the 
center  of  the  panel  except  for  boundary  condition  3  at  tha 
lower  load  levels. 

Notch  and  Cutout  Study 

The  focus  of  the  notch  analysis  involved  exaaing  the 
affect  that  different  size  cutouts  have  on  the  buckling 
behavior  of  cylindrical  panels.  Six  different  cutout  sizes 
were  studied  as  listed  in  Table  2.1.  The  2"xl”  and  l“x2" 
cutouts,  along  with  the  2"x2"  cutout  will  be  considered 
notches  for  this  thesis  because  the  area  of  the  cutout  is 
less  than  SX  of  the  total  panel  area.  Of  particular 
interest  were  the  rectangular  cutouts  with  aspect  ratios  of 
0*5  or  2.0.  The  boundary  conditions  used  for  this  part  of 
the  analysis  is  the  symmetric  boundary  condition  listed  in 
Table  2.2.  It  eay  be  recognized  that  under  the  definition 
used  in  this  thesis  to  define  a  notch,  the  resulting 
phenomena  of  collapse  tends  to  be  linear  contrary  to 
nonlinear  behavior  for  cutouts  as  will  be  shown 
subsequent! y . 
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The  procedures  developed  previously  have  been  used  to 
eodel  the  finite  element  eesh  -for  each  of  the  cutouts. 
Figure  3.33  shows  the  cutout  region’s  refined  eesh.  As  was 
noted  in  the  eesh  refinement  section  of  this  thesis*  the 
eesh  should  be  refined  equidistant  around  the  cutout.  Due 
to  the  size  and  positioning  of  the  l"x2"  and  2"xl"  notches* 
this  was  not  possible.  For  the  l"x2"  notch  the  eesh 
was  refined  2  inches  to  each  side  and  1.5  inches  to  the  top 
and  bottom  of  the  notch.  This  refinement  kept  elements  with 
a  0.5  aspect  ratio  the  desired  minimum  distance  from  the 
cutout.  Elements  with  an  aspect  ratio  of  2  were  allowed  to 
be  slightly  closer  to  the  cutout  because  their  influence  on 
the  buckling  behavior  was  not  as  critical.  In  order  to  keep 
the  elements  with  a  .5  aspect  ratio  the  einieue  distance 
from  the  opening*  the  eesh  for  the  2mxlm  notch  was  refined 
2.5  inches  in  the  circuef erential  direction.  Since  the 
2"x4"  and  4"x2“  cutouts  and  the  square  cutouts  have  even 
dimensions*  the  eesh  was  refined  evenly  around  each  cutout. 

A  suemary  of  the  grid  selection  data  for  each  of  the 
cutouts  is  presented  in  Table  3.2  The  grid  selection  data 
for  the  2"x2"  notch  with  the  19x19  mesh  can  be  found  in 
Table  3.1.  As  can  be  noted  from  Table  3.2*  the  4"x4"  cutout 
took  the  most  computer  time*  yet  had  the  lowest  collapse 
load.  This  time  function  occurred  because  there  were  over 
3000  active  degrees  of  freedom  for  this  cutout. 

The  first  parameter  examined  in  this  notch  comparison 
study  was  the  collapse  load  for  each  panel.  The  top  edge 
displacement  as  a  function  of  applied  load  is  plotted  in 


Table  3.2 


Grid  Selection  Data  for  Panels 


Size  of 
Cutout 

2"xl" 

I*x2" 

4-x4" 

2"x4“ 

4“x2" 

Grid  Size 

19x19 

17x19 

21x21 

19x21 

21x21 

Number  of 
Active  DOFs 

2761 

2437 

3007 

2903 

2917 

CPU  Time 
(sec) 

2532 

1294 

4018 

2888 

3214 

Collapse 

Load 

(lbs/in) 

292.9 

249.3 

122.0 

162.3 

180.4 

Number  of 
Load  Steps 

21 

33 

49 

32 

34 

Figures  3.34 

and  3.35. 

The 

panels  with 

the  2-xl“ 

and  l"x2 

notches  have  higher  collapse  loads  because  there  is  less 
eaterial  cut  f roe  the  panel.  The  2"xl"  notch  produces  the 
highest  collapse  load  because  its  seal lest  direction  is 
perpendicular  to  the  applied  load.  This  is  a  stiff er  panel 
than  the  l"x2"  notch  panel  because  the  l"x2"  notch  has 
higher  stress  concentrations  at  its  corners.  The  higher 
stress  concentrations  can  be  correlated  to  higher  strains 
Mhich  results  in  greater  displacements.  The  4“x4"  cutout 
has  the  lowest  collapse  load  because  it  has  the  most 
material  removed.  The  4"x2"  cutout,  which  has  an  aspect 
ratio  of  2V  has  a  higher  collapse  load  than  the  2“x4a 
cutout.  This  is  due  again  to  the  high  stress  concentrations 
at  the  corners  of  the  cutout  when  considering  the  0.5  aspect 
ratio.  These  figures  show  that  notches  with  an  aspect 
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Figure  3.33.  Load  Displacement  Curves  for  Panels  with 
2"k4",  4"k2",  and  4"x4"  Cutouts 
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ratio  of  2  can  withstand  greater  loads  than  cutouts  of  the 
saee  area  with  aspect  ratios  of  either  0.5  or  1. 
Furthermore,  these  figures  show  that  those  cutouts  defined 
as  notches  resist  load  linearly  relative  to  displacement. 

As  mentioned  previously,  the  area  around  cutouts  with 
0.5  aspect  ratios  deflect  more  than  notches  of  the  same  area 
with  aspect  ratios  of  2  at  equivalent  loads.  This  can  be 
seen  in  Table  3.3  which  lists  the  maximum  and  minimum  radial 
displacements  corresponding  to  applied  loads.  For  the 
2a>xl",  l"x2"  and  2"x2"  notches,  the  point  of  minimum  radial 
displacement  is  basically  the  same.  It  occurs  at  the  center 
of  the  upper  edge  of  the  cutout  in  all  cases.  The  point  of 
maximum  radial  displacement  shifts  from  the  upper  right  edge 
of  the  cutout  at  lower  load  levels  to  the  upper  left  quarter 
of  the  panel  at  collapse  load.  This  same  pattern  can  be 
seen  for  the  4"x2'*,  2"x4"  and  4"x4"  cutouts.  Minimum  radial 
displacement  occurs  near  the  center  of  the  top  edge  of  the 
cutout.  Maximum  radial  displacement  at  collapse  load  is  at 
the  upper  left  quarter  of  the  panel.  The  2"xl"  notch 
produces  the  least  amount  of  radial  displacement  at  collapse 
load.  The  panel  with  the  4"x4"  cutout  has  the  greatest 
amount  of  surface  area  removed  which  results  in  a  more 
flexible  structure.  Therefore,  this  panel  has  the  largest 
radial  displacement  at  collapse  load  of  all  the  panels 
tested . 

To  analyze  the  displacement  pattern  over  the  entire 
panel,  contour  plots  at  collapse  load  are  presented  in 
Figures  3.36  through  3.41.  For  the  2"xl",  l-x2"  and  2"x2" 
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notch  plots*  there  is  a  definite  inward  deflected  trough  of 
radial  di spl acaaant  which  axtsnds  froa  tha  uppsr  left  adga 
of  tha  cutout  to  tha  lowar  right  adga.  For  tha  2“xl“  notch* 
tha  trough  is  alaost  diagonal  froa  cornar  to  corner  whereas 
tha  trough  for  tha  l"x2“  and  2Mx2H  begins  to  shift  aora 
towards  tha  center  of  the  panel.  Two  bulges  of  positive 
radial  di spl acaaont  occur  along  either  side  of  tha  trough 
with  the  aaxiaua  di spl acaaant  being  on  tha  upper  laft 
quarter  of  the  panel.  Refer  to  Table  3.3  for  tha  coordinate 
locations  of  aaxiaua  and  ainiaua  radial  di spl acaaant.  For 
tha  4“x2",  2"x4"  and  4"x4"  cutout  contours  (Figures  3.39 
through  3.41)*  this  trough  does  not  extend  diagonally  but  is 
basically  straight  up  and  down  tha  panel.  The  di spl acaaant 
pattern  of  tha  panels  with  the  4“x2"  and  2“x4“  cutouts 
reaain  alaost  syaaetric  at  collapse  load  and  tha  panel  with 
the  4"x4"  cutout  raaains  entirely  syaaetric  at  collapse.  As 
in  the  other  contour  plots  the  width  of  the  trough 
corresponds  to  the  width  of  the  cutout.  Janisse  C231 
discusses  the  fact  that  the  linear  eigenvector  for  a  2"x2" 
cutout  is  siailar  to  the  collapse  contour  pattern*  thus  the 
saall  cutout  reacts  linearly.  This  is  not  true  with  the 
larger  cutouts.  This  verifies  what  was  stated  previously 
for  the  definition  of  a  notch  in  this  thesis. 

To  further  analyze  the  effects  of  cutout  size  on  the 
collapse  behavior  of  the  panel*  radial  displacement  and 
aoaent  about  the  y  axis  were  analyzed.  Radial  displacement 
as  a  function  of  distance  from  the  left  edge  of  the  panel  is 
plotted  for  the  lHx2"  notch  in  Figure  3.42.  Again  the 


wa<|Tt  ■  .02547  inchn 

w  .  -  -.05377  inchn 

•in 

Figure  3.37.  Radial  Displacement  Contour  Plot  for  l"x2" 
Notch  at  Collapse  Load  of  249.3  lbs/in 
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Not*:  Positive  radial  displacaaant  ia  outward 


Figura 


WJBax  ■  .02143  inchas 
••  4  -  -.04391  inchas 

.38.  Radial  Displacement  Contour  Plot  for  2"x2 
Notch  at  Collapse  Load  of  207.3  lbs/in 


Mots:  Positive  radial  displacement  is  outward. 

w(tiw  *  .07917  inches 

w  ,  -  -.04535  inches 

mi  n 

Figure  3.39.  Radial  Displacement  Contour  Plot  for  4"x2 


Cutout  at  Collapse  Load  of  180.4  Ibs/in 


Not*:  Positive  radial  displace want  is  outward. 

wwMT  ■  .#9372  inches 
■  -.08487  inches 

Figure  3.41.  -  Radial  Displaceeent  Contour  Plot  for  4*k4" 
Cutout  at  Collapse  Load  of  122.0  lbs/in. 
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trough  begins  to  appear  with  increasing  load  and  the  aini 


radial  di spl aceeent  shifts  froe  about  4"  f roe  the  left  edge 
(along  the  line  3"  froe  the  top  of  the  panel)  to  about  7" 
froe  the  left  edge  (along  the  line  9"  froe  the  top  of  the 
panel).  The  aoasnt  about  the  y  axis  as  a  function  of 
distance  froe  the  left  edge  is  plotted  in  Figures  3.43 
through  3.45  for  the  laax2n  notch.  The  point  of  ainiaua 
along  a  line  4.25"  froe  the  top  of  the  panel  corresponds  to 
the  saee  location  as  the  point  of  ainiaua  radial 
di spl aceeent  in  Figure  3.42.  The  eoeents  are  greater  in 
aagnitude  along  the  area  where  the  trough  fores  in  the  panel 
than  near  the  right  edge  of  the  panel.  Figure  3.44  is  a 
plot  of  Hx  along  a  line  5.75"  froe  the  top  of  the  panel 
•diich  passes  through  the  cutout.  At  collapse  load,  there  is 
a  Juap  froe  negative  aoaent  on  the  left  edge  of  the  notch  to 
•  larger  positive  aoaent  on  the  right  edge  of  the  notch. 
Along  the  line  7.75"  froe  the  top  of  the  panel  (Figure 
3.45),  the  eaxiaua  negative  eoeents  occur  in  the  saee  area 
that  the  trough  fores.  At  lower  load  levels,  the  values  for 

*r»  syeeetric  along  each  of  the  lines  across  the  panel. 
At  collapse  load,  the  eoeents  are  not  syeeetric  and  the 
largest  negative  eoeents  occur  along  the  radial  displaceaent 
trough. 

The  radial  displaceaent  and  aoaent  in  the  x  direction 
for  the  2"xl"  notch  are  plotted  in  Figures  3.46  through 
3.49.  The  radial  displaceaent  trough  seen  in  Figure  3.46 


extends  froe  approxiaately  5"  froe  the  left  edge  (above  the 
cutout)  to  7”  froa  the  left  edge  (along  the  line  below  the 


•V 


Figure  3.42.  Radial  Displacement  Profile  at  Collapi 
Load  for  lMx2"  Notch 
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.43.  Mk  along  a  Lina  7. 75"  from  Top  Edg 


of  Panal  for  l”x2"  Notch 


T  (INCHES) 

0.01  0.02  0.03  0.04  0.05 
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2.0 

DISTANCE 
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ALONG  LINE  3'  FROM  TOP  EDGE 
ALONG  LINE  6'  FROM  TOP  EDGE 
ALONG  LINE  9‘  FROM  TOP  EDGE 


Figure  3.46.  Radial  Displacement  Profile  at  Collapse 
Load  for  2"xl"  Notch 


cutout).  The  notch  is  slightly  skewed  because  ths  left 
•dgs  of  ths  notch  displaces  inward  sore  than  the  right  edge 
of  ths  notch.  Positive  radial  displaceeents  correspond  to 
positive  eoeents  in  Figures  3.47  through  3.49.  Hinieum 
noeents  correspond  to  the  locations  of  einiaua  radial 
displaceeent  (the  trough)  in  Figure  3.46.  Along  the  line 
4.25“  f roe  the  top  of  the  panel,  the  eoaent  is  snail  near 
the  right  edge  of  the  panel  and  the  einiaua  value  occurs  in 
the  radial  displaceeent  trough.  The  eoeents  are  priaarily 
positive  along  the  line  that  passes  through  the  notch 
(Figure  3.48).  On  the  left  edge,  the  aoment  goes  negative 
at,  collapse  load,  then  eakes  a  large  juap  to  a  positive 
eoaent  on  the  right  edge  of  the  notch.  This  is  due  once 
again  to  the  cutout  being  skewed  as  load  is  applied.  The 
values  for  are  very  seall  along  the  line  7.75"  fore  the 
top  (Figure  3.49)  except  for  a  large  negative  eoaent  in  the 
area  of  the  trough. 

The  radial  displaceeents  and  eoaent  resultants  for  the 
2"x2“  notch  are  presented  in  Figures  3.50  through  3.53.  The 
positive  radial  displacements  in  Figure  3.50  correspond  to 
positive  values  of  Along  the  line  4.25“  from  the  top  of 
the  panel  (Figure  3.51),  the  moments  near  the  right  edge  of 
the  panel  are  less  than  near  the  left  edge.  This  again  can 
be  attributed  to  the  buckling  shape  of  the  panel  with  the 
trough  forming  from  the  upper  left  section  of  the  panel  to 
the  lower  right  section.  The  moments  directly  above  the 
notch  are  negative.  Along  the  line  that  passes  through  the 
notch  Figure  3.52,  the  moment  again  jumps  in  magnitude  from 
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Figure  3.47.  along  a  Lin*  4.25"  from  Top  Edg* 


of  Panel  for  2"xl“  Notch 
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Figure  3.50.  Radial  Displacement  Profile  at  Collapse 


Load  for  2"x2”  Notch 
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Figure  3. Si.  along  a  Line  4.25“  from  Top  Edg 

of  Panel  for  2"k2"  Notch 
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Figura  3. 92.  along  a  Lina  5.75"  from  Top  Edg 

of  Panal  for  2”x2N  Notch 
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the  left  edge  to  the  right  edge  of  the  notch.  The  largvst 
■onnti  in  the  panel  at  col  laps*  load  ara  at  tha  right  adga 
of  tha  cutout.  Thara  is  a  larga  nagativa  aoaant  in  tha 
trough  along  tha  7.75"  lina  (Figura  3.53)  which  corrasponds 
to  tha  larga  nagativa  aoaant  in  tha  trough  along  tha  uppar 
lina. 

For  tha  2"x4"  cutout ,  tha  radial  displacaaant  trough  is 
not  diagonal  across  tha  panal  as  with  tha  saallar  notchas. 
Tha  trough  is  also  considerably  wider,  corresponding  to  tha 
width  of  tha  cutout.  Figura  3.54  depicts  tha  displacaaant 
pattern  for  tha  2"x4"  cutout.  Tha  point  of  ainiaua  radial 
displacaaant  shifts  to  the  right  by  only  ana  inch  froa  top 
to  bottoa  of  tha  panal.  For  low  load  levels,  is 
syaaetric  about  tha  canter  of  tha  panal  as  can  be  seen  in 
Figura  3.55.  At  collapse  load,  tha  aoaant  still  retains  its 
syaaatry  along  the  uppar  lina.  The  ainiaua  value  for  tha 
aoaant  once  again  occurs  in  tha  slight  radial  displacaaant 
trough.  along  the  line  at  tha  canter  of  tha  panal, 
Figura  3.56,  is  greater  at  tha  left  side  of  tha  cutout. 
This  is  apposite  to  what  happens  with  tha  smaller  notchas. 
This  affect  can  be  attributed  to  the  fact  that  there  is  not 
a  deep  radial  displacaaant  trough  extending  diagonally 
across  tha  panal  contrary  to  tha  smaller  notchas.  Tha 
moments  along  tha  lina  near  tha  bottom  edge  (Figura  3.57) 
ara  less  than  tha  moments  at  other  locations  and  retain 
thair  symmetry  throughout  tha  load  application  process. 
This  is  due  to  tha  bottom  line  being  a  large  distance  from 
the  cutout  where  the  affect  of  the  large  displacements 
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Figure  3.54.  Radial  Displacement  Profile  at  Collapse 


Load  for  2"k4”  Cutout 


MX  (IN-LBS/IN) 

-4.0  -3.0  -2.0  -1.0  0.0  1.0  2.0  3.0  4.0 


•round  the  cutout  arc  not  felt 


Also,  this  lins  is  ncari 


ths  clamped  bottom  edge  which  restrains  displacements  and 
therefore  results  in  smaller  moments. 

Once  again  the  radial  displacement  profile  (Figure 
3.58)  shows  only  a  slight  shift  in  the  location  of  the 
displacement  trough  from  top  to  bottoe  of  the  panel  for  the 
4"x2"  notch.  There  are  large  symmetric  peaks  in 
displacement  on  either  edge  of  the  cutout.  This  is  due  to 
the  small  length  of  the  notch  in  the  circumferential 
direction  which  causes  the  notch  not  to  skew  with  applied 
load.  The  moments  along  the  3.75"  line  are  much  greater 
than  the  moments  for  the  smaller  notches  as  can  be  seen  in 
Figure  3.59.  The  moments  are  symmetric  about  the  center  of 
the  panel  for  all  load  levels.  This  is  because  the  trough 
is  not  diagonal,  but  rather  follows  the  center  line  of  the 
panel.  along  the  center  of  the  panel  (Figure  3.60)  is 

syeuaetric  about  the  notch,  again  indicating  that  the  cutout 
does  not  skew.  Along  the  line  near  the  bottom  of  the  panel. 
Figure  3.61,  the  moments  are  greater  than  for  the  smaller 
notches  and  are  symmetric. 

The  results  for  the  4"x4"  cutout  follow  the  same 
pattern  as  for  the  4"x2“  cutout.  Figure  3.62  shows  the 
radial  displacement  trough  does  not  shift  significantly  to 
the  right.  Maximum  radial  displacement  occurs  along  a  line 
near  the  center  of  the  cutout.  For  the  panel  with  the  4"x4" 
cutout,  the  radial  displacement  plot  shows  that  there  are 
large  displacements,  and  consequently,  large  rotations  near 
the  cutout  at  collapse  load.  The  maximum  displacements  are 


MX  { IN-LBS/IN) 

-8.0  -6.0  -4.0  -2.0  0.0  2.0  4.0  6.0  8.0 


DISTANCE  PROM  LEFT  EDGE  (INCHES) 


LEGEND 

□  -  60.6  LB  LOAD 
a  -  128.5  LB  LOAD 
o -  178.0  LB  LOAD 


Figure  3.58.  Radial  Displacement  at  Collapse  Load 
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Figure  3.S9.  along  a  Line  3. 75"  from  Top  Edge 


of  Panel  for  4nx2M  Cutout 
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Figure  3.60.  Along  a  Line  3.75"  from  Top  Edgi 

of  Panel  for  4"x2"  Cutout 
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Figure  3.62.  Radial  Displacement  Profile  at  Collapi 
Load  for  4"x4"  Cutout 
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over  twice  the  thickness  of  the  panel.  Since  a  large  aeount 


of  the  surface  area  of  the  panel  is  reeoved,  the  edges  of 
the  cutout  act  like  free  edges  and  displace  large  aeounts. 
This  panel  and  the  panel  with  the  4Mx2*  cutout  are  the  first 
cutouts  exaained  in  this  thesis  that  truly  eake  use  of  the 
large  displaceeent  theory  portion  of  STAQSC— 1.  The 
rotations  are  even  larger  than  theory  Indicates  they  should 
be.  The  large  rotation  characteristics  in  the  vicinity  of 
the  cutout  are  also  due  to  the  free  edge  effect.  The  ooeent 
along  each  of  the  lines  (Figures  3.63  through  3.65)  are  all 
syeeetric  about  the  center  of  the  panel.  The 
characteristics  of  the  eoeents  are  different  f roe  those  for 
the  notches.  The  eoeents  for  these  large  cutouts  are 
syeeetric  about  the  cutout  at  collapse  load.  Again,  Janisse 
C233  pointed  out  that  large  cutouts  (4"x4">  displace 
differently  when  coeparing  the  nonlinear  analysis  with  the 
linear  results.  The  eigenvector  at  buckling  contains  a 
trough  sieilar  to  the  seall  notch  cutouts  previously 
considered.  The  collapse  displaceeent  patterns  for  the 
large  cutouts  do  not  follow  this  pattern,  but  are  relatively 
syeeetric  which  is  a  characteristic  of  the  ieportance  of 
bending  in  the  resistance  for  the  structure  against  bending. 
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IV. 


USIONS 


B*«td  on  th«  analysis  prsssntsd  in  this  thesis,  tha 
following  conclusions  can  ba  aada: 

1.  Individual  alaaants  in  a  finita  alaaant  aash  should 
praferably  have  equal  dimensions  in  tha  x  and  y  directions. 

2.  In  a  refined  mesh  around  the  area  of  a  cutout, 
rectangular  QUAF  411  elements  influence  the  moments  in  the 
direction  of  their  longer  dimension  to  a  greater  degree  than 
square  elements.  If  elements  near  the  cutout  have  an  aspect 
ratio  of  0.5f  the  collapse  load  estimation  is  not  improved 
as  one  would  expect  over  the  results  for  a  coarse  mesh. 

3.  In  a  mesh  refinement  process  from  one  inch  square 
elements  to  half-inch  square  elements,  the  half-inch  square 
elements  should  be  a  uniform  distance  around  the  cutout  for 
a  STA6SC-1  nonlinear  collapse  analysis.  Elements  with  an 
aspect  ratio  of  0.S  should  be  at  least  two  inches  from  the 
cutout  and  elements  with  an  aspect  ratio  of  2  should  be  at 
least  one  and  a  half  inches  from  the  cutout. 

4.  Restraining  the  v  displacement  along  the  upper  left 
or  upper  right  half  of  the  cylindrical  panel  does  not 
significantly  improve  the  results  over  the  results  obtained 
with  symmetric  boundary  conditions  when  compared  to 
experimental  data. 

5.  Experimental  results  considering  the  given  panel 
and  supports,  are  better  approximated  by  an  analytic 
boundary  condition  in  which  the  u  and  v  displacements  are 
held  fixed  along  the  left  edge  of  the  panel  between  1  and  6 
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inches  -from  the  top  edge  of  the  panel.  This  boundary 
condition  appears  to  eodel  the  change  in  boundary  conditions 
that  occurs  experieentally  with  applied  load. 

6.  Surface  imperfections  in  the  test  panel  result  in  a 
lower  collapse  load  and  different  radial  displacement 
patterns  than  those  modelled  analytically. 

7.  The  restrained  boundary  condition  introduces 
greater  stiffnesses  into  the  panel  which  result  in  a  higher 
collapse  load. 

8.  It  appears  to  be  possible,  with  the  use  of  the 
STAGSC-1  computer  program,  to  actually  model  an  experimental 
set  of  boundary  conditions.  The  ideal  boundary  condition 
was  not  reached  in  this  thesis,  but  with  additional  computer 
time,  several  iterations  could  be  carried  out  and  the  ideal 
boundary  condition  attained. 

9.  For  a  complete  numerical  analysis  of  a  panel  tested 
experimentally,  one  .  oust  take  into  account  surface 
eccentricities  as  well  as  varying  boundary  conditions. 

10.  Panels  with  notches  having  an  aspect  ratio  of  2 
have  higher  collapse  loads  than  panels  with  notches  of  the 
same  area  but  with  aspect  ratios  of  0.S.  This  is  due  to  the 
higher  stress  concentrations  around  the  corners  of  the 
cutouts  with  0.5  aspect  ratios.  The  higher  stresses  result 
in  larger  displacements  and  therefore  lower  collapse  loads. 
Panels  with  cutouts  with  an  aspect  ratio  of  0.5  can  carry  up 
to  17X  greater  load  than  panels  with  the  same  area  cutout 
having  an  aspect  ratio  of  2.0. 

12.  The  2"xi"  notch  is  associated  with  the  highest 
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collapse  load  but  uullMt  radial  displacements.  This  is 
because  tha  notch  is  a  saall  par cant aga  of  panal  araa  and 

theref ore,  tha  panal  ratains  auch  of  tha  structural 

intagrity  of  a  panal  without  a  cutout. 

13.  Tha  aoaants  in  tha  x  diraction  ara  syaaatric  about 

all  cutouts  at  low  load  lavals.  At  collapsa  load,  tha 
aoaants  for  the  4Mx2>>  and  4"x4"  retain  thair  syaaatry 

indicating  that  these  cutouts  do  not  skew  with  applied  load 
as  is  the  case  with  tha  other  cutouts. 

14.  Considering  the  author’s  definition  of  a  notch,  (a 

cutout  of  lass  than  51  of  the  panel  area)  panels  with 
notches  behave  differently  froa  a  collapse  displacement 
function  point  of  view  than  panels  with  cutouts.  The  radial 
displaceaent  pattern  of  a  panel  with  a  notch  foras  a 

diagonal  trough  whereas  panels  with  cutouts  have  a  vertical 
displacement  trough  down  the  center  of  the  panel.  Thus,  the 
notch  does  not  produce  tha  typical  nonlinear  effect 
associated  with  cutouts. 
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